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We consider application of domain wall fermions to quarks with relatively heavy masses, aiming at precision 
calculations of charmed meson properties. Preliminary results for a few basic quantities are presented. 



1. INTRODUCTION 



2. SIMULATION PARAMETERS 



After the remarkable accomplishment of the 
B factories, theoretical uncertainty remaining in 
the decay amplitudes of heavy-light mesons is an 
important obstacle to the stringent test of the 
Standard Model. Lattice QCD offers a promis- 
ing method to reduce this uncertainty in a sys- 
tematic way and numerous attempts have been 
made to achieve this purpose so far |1I2| . While 
the accuracy we can attain in simulating b quark 
is still limited, mainly due to difficulty in tak- 
ing a continuum limit, present computational re- 
sources should make it possible to perform preci- 
sion calculations of charmed hadron properties. It 
is expected that the charm factories will allow us 
to calibrate lattice calculations of the heavy-light 
system with high precision, and the experiences 
gained in c-physics will be advantageous to the 
&-physics simulations. 

Domain wall fermions (DWF) |3I4I5| are ex- 
pected to have only a negligible 0{a) error as well 
as better control over chiral extrapolations. Non- 
perturbative renormalization, which is essential 
to high precision calculations, also seems to work 
remarkably well j6j. Simulating light quarks us- 
ing DWF benefits from these advantages and has 
been successfully used in the light hadron system 
|7IXI9| . It is, then, natural to consider its ap- 
plication to the massive quarks. As a first step 
toward the precision calculations, we have started 
the simulation of D mesons using DWF for both 
heavy and light quarks. 



The first exploratory study of massive DWF 
|1U| observed that with an improved gauge ac- 
tion at weaker gauge coupling DWF do not show 
clear failure for anig'^OA. Based on this ob- 
servation, we take the DBW2 gauge action [TT] 
with /3=1.22 at which the charm quark is realized 
around amq ~0.4. The simulation is carried out 
using quenched lattices with 24^ x 48. We set the 
extension of the fifth dimension and the domain 
wall height to Lg^lO and a7\f5=1.65, respec- 
tively. The light and heavy quark masses range 
in 0.008< amiq <0.040 and 0.1< am,,, <0.5, re- 
spectively. To avoid the finite size effect in the 
time direction, all the meson two-point functions 
are calculated under both the periodic and anti- 
periodic boundary conditions and then averaged. 
The gauge field is fixed to the Coulomb gauge 
in this work, and a~^=2.86(9) GeV is obtained 
using nip— 770 MeV. Data presented below is ob- 
tained with 20 configurations. 

The magnitude of the explicit chiral symme- 
try breaking can be probed by the residual mass, 
TTircs, defined by 

where the explicit forms of the operators are 
found in Ref. j7l9j . and we obtain arrives = 
9.72 X 10^^, which roughly corresponds to mrcs ~ 
0.3 MeV. 

3. PRELIMINARY RESULTS 



*We thank RIKEN, Brookhaven National Laboratory and 
the U.S. Department of Energy for providing the facihties 
essential for the completion of this work. 



To see if the large mass parameters cause any 
unexpected behavior, we first examine the effec- 
tive mass plots. The plots shown in Figure Q are 
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Figure 1. Effective mass plots of the heavy-light 
pseudo-scalar mesons. The fit ranges and the re- 
sults are indicated by horizontal lines. 
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Figure 2. Light quark mass dependences of the 
heavy-light pseudo-scalar (circles) and vector (tri- 
angles) meson masses. 



obtained from the {A4P^) correlation functions 
with wall source and point sink, where A4 and P 
are heavy- light bilinears. The light quark mass is 
fixed to aTO/(j= 0.040 while the heavy quark mass 
varies from 0.1 (bottom) to 0.5 (top) by 0.1. The 
data with aruhq^OA and 0.5 show a gradual in- 
crease at t^3Q and it looks getting more conspic- 
uous as aruhq gets larger, otherwise no odd be- 
havior is seen for all heavy quark masses. We re- 
peated the same calculation for one configuration 
with more stringent stopping condition and found 
that in the heaviest case {amhq—0.5) the corre- 
lation function with the smaller stopping condi- 
tions differs by only -fO.05% at t=20 while the 
difference reaches to -t-10% around t ~30. There- 
fore it seems that the increase found in Figure ^ 
is caused by the loose stopping condition. Since 
all the correlation functions for t <20 show no 
odd behavior and are seen to have a reasonable 
plateau, we take this range to perform a conser- 
vative analysis. The errors given below are sta- 
tistical only, unless otherwise specified. 

Figure [5] shows the chiral extrapolation of the 
heavy-light pseudo-scalar (circles) and vector (tri- 
angles) meson masses with amiiq=OA. Both 
the linear and quadratic fits describe the data 
very well, and the difference due to the dif- 
ferent fits does not exceed 0.5 % in the chi- 



ral limit for the both mesons. We choose the 
linear fit to obtain the heavy-light masses at 
a'mud=~o,m,-cs and am^^O. 032(2) for later use, 
where amg, corresponding to the strange quark 
mass, is determined using rriK/mp as a rep- 
resentative among other possible inputs. To 
find amcharmi the pseudo-scalar masses are in- 
terpolated to the physical value of tod using a 
quadratic function of heavy quark mass, and we 
obtain amcharm=0. 380(4). The mass splitting of 
niD, — niD is, then, obtained as 

mo, -mo^ 114(6) MeV. (2) 

While this value is about 15% larger than the 
experimental value (99MeV), quenching may ac- 
count for much of the discrepancy. 

The lattice calculation of the IS hyperfine 
splitting has been a long standing problem be- 
cause the lattice results are significantly smaller 
than the experimental values, independent of the 
heavy quark action adopted. In this work we ob- 
tain 

mD--mD = 69(6) MeV, (3) 
mo; - mo, = 70(4) MeV. (4) 

Comparing to the experimental results of mjy — 
r7i£i = 142 MeV and mo* -r7iD,=144 MeV, we see 
that this problem persists when using DWF. 
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Figure 3. The chiral extrapolation of the pseudo- 
scalar decay constant at amhq=0.4:. 



The leptonic decay constants are obtained in a 
standard way from two-point functions with wall 
source and wall sink. Such correlation functions 
are still rather noisy, and the resulting decay con- 
stants depend on the fit ranges by as much as 
10%. In the following analysis, we thus adopt 
a set of fit ranges and take 10% uncertainty as a 
systematic error. Figure|31shows the chiral behav- 
ior of the decay constant at amhq=OA. Since the 
data with am;g=0.008 looks unnaturally high, we 
performed linear extrapolations with and with- 
out this data point. The difference between the 
two extrapolations are taken as a systematic er- 
ror again. After interpolation to amhq—amcha.vm, 
we obtain 



/•latt 

Jd 



j'la 

Jd 



latt 



200(9)(l2i) MeV, 
216(7)(22) MeV, 



(5) 
(6) 



where the systematic errors are summed up in 
quadrature. At the moment we quote the lattice 
bare values because we do not know how signifi- 
cant the 0(a"TO"ijarm) error is. 

The SU(3) flavor breaking ratio of the decay 
constant is one of the quantities for which lattice 
QCD can attain a high precision because most of 
uncertainties are expected to cancel in the ratio. 
We can anticipate that the on-going charm ex- 
periments make it possible to measure this ratio 
at a few % accuracy, and so this is a good quan- 



tity for which to compare the lattice result with 
the experimental one although the uncertainties 
in nig , the chiral extrapolation and quenching ef- 
fects must first be under control. Our current 
result is 



fDjfD = 1.08(3)(tg) 



4. SUMMARY 



(7) 



We have started a exploratory study of mas- 
sive DWF aiming at precision lattice calculations 
of the D meson system. The preliminary results 
show no evidence of the DWF mechanism failing 
at large masses, while it may be better to im- 
prove how to impose a stopping condition. How- 
ever, for precision calculations, we need to un- 
derstand scaling violations, especially how seri- 
ous the 0(a"TO"jj3^j.jjj) error is. We also need more 
statistics to have better control over chiral ex- 
trapolations. 
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